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A dissipative mechanism is presented, which emerges in generic interacting 
quantum field systems and which leads to robust warm inflation. An explicit 
example is considered, where using typical parameter values, it is shown that 
considerable radiation can be produced during inflation. The extension of our 
CN ' results to expanding spacetime also is discussed. 
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I. INTRODUCTION 



Inflationary dynamics inherently is a multifield problem, since the vacuum energy that 
drives inflation eventually must convert to radiation, which generally is comprised of a 
^ ■ variety of particle species. Phenomenologically it has been shown that the inflation and ra- 
diation production phases can be two well separated periods in scenarios generically termed 
supercooled (or isentropic) inflation (for a review see [1]), or radiation production can occur 
concurrently with inflationary expansion in scenarios generically termed warm (or nonisen- 
tropic) inflation [2]. Warm inflation is a broader picture, since the extent of radiation 
production during inflation is variable, so that supercooled inflation emerges as the limiting 
case of zero radiation production. 

Although by now considerable work has demonstrated its phenomenological significance 
[3], one key barrier to the warm inflation picture has been establishing plausibility of its 
dynamics from first principles quantum field theory. To some extent this point has been 
overemphasized for warm inflation, since in similar respects particle production during the 
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far out-of-equilibrium reheating phase of supercooled inflation is not well understood, thus 
leaving incompleteness also to this picture. However for supercooled inflation, since particle 
production is assumed not to affect large scale structure formation during inflation, thus the 
main observational predictions, these shortcomings are cast aside as secondary concerns. 
Nevertheless, without a solution here, this picture is unproven. On the other hand, the 
warm inflation picture makes no a priory assumption that particle production does not 
affect large scale structure formation. As such, the particle production problem appears 
more acute here. More basically a proper understanding of particle production should mean 
that theory itself can decide which or to what extent either of these two pictures is valid. 
Undoubtedly, no theory based on inflationary expansion will ever emerge, until particle 
production in quantum field theory is adequately understood. 

This is a major problem, which must be tackled in steps. Fair enough is to attempt to 
see how well either picture of inflation can be understood from first principles and en route 
hope a clearer general picture eventually will emerge. For warm inflation, there is greater 
possibility to understand particle production, and eventually reach closure at a theoretical 
level about the viability of this picture as a description of the early universe. The reason 
is that recall in this picture the scalar inflaton field is required to have a slow, overdamped 
motion. As such, adiabatic methods of quantum field theory are applicable here, and these 
are the only methods for which dissipation can be unarguably analyzed. 

The road toward a first principles warm inflation picture primarily has been hindered 
by basic gaps in the understanding of dissipative quantum field theory, which during the 
course of developing warm inflation are being filled [4-9]. The first attempt to understand 
warm inflation dynamics utilized finite temperature dissipative quantum field theory, since 
some formalism already existed here [10-14]. Based on this work [4], statements of a general 
sort have been made about the impossibility of warm inflation dynamics [6]. However, these 
criticisms failed to recognize that the key problems were specific to the restrictive constraints 
of the high-T approximation and were not reflexive of warm inflation in general. 

Intrinsically, warm inflation is an out-of-equilibrium problem, in that it is not tied to any 
specific equilibrium statistical state, but rather simply requires radiation production con- 
current with the overdamped relaxation of a global order parameter. Although the actual 
statistical state during warm inflation may not be very far from an equilibrium state, at 
present the problem is simply technical limitations in describing the scope of such states. 
Furthermore, as has been noted [2,7], very little radiation production during inflation, at 
the scale of tens of orders of magnitude below the vacuum energy density, is already suffi- 
cient to affect large scale structure formation and create an adequately high post-inflation 
temperature. 

With these thoughts in mind, in [7] a simple attempt was made to circumvent the 
specific constraints of the high-temperature formalism, by examining dissipation at zero- 
temperature. The point there was to investigate a suggestions learned from our high- 
temperature analysis, that alleviation of the constraints specific to the high-T approximation 
would adequately allow realizing robust radiation production during warm inflation. The 
main purpose of [7] was to develop the necessary formalism, but in addition one suggestive 
mechanism was identified that could realize this point, which involved a scalar $ field (whose 
zero mode can be associated, e.g., with the inflaton) exciting heavy x-bosons which then 
decay into lighter -?/>-fermions. This letter reports a detailed investigation of this process 
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and demonstrates that it is a robust mechanism for warm inflation. For this, in Sec. II a 
linear response derivation will be presented, which in the adiabatic regime and at leading 
order is equivalent to the closed time Lagrangian formalism, but is simpler and physically 
more transparent. Then in Sec. Ill an alternative derivation is presented, using canonical 
methods. From this approach, the origin of particle production and energy balance for this 
mechanism will be clarified. Next, Sec. IV gives a physical picture to the mechanism and 
supplies an explicit numerical example to demonstrate the extent of radiation production 
it yields during inflation. Sec. V discusses the extension of the calculation to expanding 
spacetime. Finally the conclusions are given in Sec. VI. 



II. A MODEL FOR ROBUST RADIATION PRODUCTION 

We consider a multi- field model, first studied in [7], of a scalar field $ interacting with 
a set of scalar fields Xj, 3 — N x , which in turn interact with fermion fields ipk, k — 

1, . . . , N^, with Lagrangian density 



\ 2 ^^Ox ' 2 



4] A 3 -i 



fc=l 



0-m OlPk -J2 h kjXj 
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(1) 



The regime of interest for warm inflation, that is studied here is m Xj > 2m 1 p k > m^, where 
these are the renormalized and, if relevant, background field dependent masses. 

By decomposing $ in terms of a homogeneous classical part, <p(t), and its fluctuations 
0, the effective equation of motion (EOM) for (p emerges as 



A 



A 



A 



<p(t) + mlMt) + -ip\t) + -<^(t)(0 2 ) + -(0 



N x 
3=1 



(2) 



We will use a linear response theory approach in which the field averages in Eq. (2) are 
expressed in terms of the respective field propagators G,p(x,x') and G Xj (x,x'). Also in the 
following, we derive the tp effective EOM from an adiabatic approximation. This approxima- 
tion requires that all macroscopic motion is slow relative to the characteristic scales of the 
microscopic dynamics. In our model the time scale for microscopic dynamics is represented 
through the (inverse of the) particle decay widths T^, T x and for macroscopic dynamics is 
contained in <p(t), with the basic consistency condition [4] 

vV^«r^,r x . (3) 

Turning to the derivation, consider first (Xj)- This expectation value can be expressed 
in terms of the coincidence limit of the (causal) two-point Green's function for the Xj field, 
G x + (x, x'). Recall that this Green's function is the (1, l)-component of the real time matrix 
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of full propagators, all of which satisfy the appropriate Schwinger-Dyson equations (see, e.g., 
[4,7] for additional details) 



□ + m\ . + g?<p 2 (t)\ G Xj (x, x' 



+ J d*zE Xj (x, z)G Xj (z, x') = i5(x, x') , (4) 

where £ Xj . is the Xj field self-energy. The field frequencies appearing in these propagators 
depend on the background field configuration (p(t). This field is decomposed as (p(t) = 
y?o + $¥{t), where <po is a constant (the value of the field at say the initial time t = to) and 
5<p(t) is treated perturbatively. This is just a linear response theory approach to calculating 
the averages of the fields appearing in Eq. (2). Following this procedure, we have that (x 2 ) 
can be written to lowest order as 



-i f dt' Q l [p\t')-^ <[^(x,t), Xj 2 (x,t')]> , (5) 

J -co Z 



where (. . .)o means the correlation function evaluated at the initial time. The (p 2 dependence 
in Eq. (5) emerges from expanding the two point function with respect to the 8tp dependent 
terms. Formally this can be done by treating Sip dependent terms in the shifted potential 
as interaction vertices. This implies adding an interacting vertex quadratic in the Xj field, 
with Feynman rule —ig 2 /2 [(f 2 (t) — (fl], and is used in calculating the leading order one- loop 
bubble diagram that gives the two-point function. This method was first implemented to 
study dissipation in [12,13] and more recently in [7]. This is also analogous to the functional 
Schwinger closed time path formalism used in [10,4]. Using translational invariance we can 
now write ([x|(x, t), x|(x, f)]), appearing in Eq. (5), in terms of the two-point Green's 
function for the Xj field, G++(x, x'), as 

(bd(x, t), X 2 (x, t')}) = 2i Im(T X Kx, t) X 2 A^ f )> 



4 * / T^M^J+fa, t - t')f t>t , , (6) 



v 3 

Al I (2tt) 3 "" l "^ 

where G+ + (q, t — t') is given by (see e.g. [7] for the explicit expressions for both the scalar 
and fermion field propagators) G++(q, t-t') = G>. (q, t-t')9(t- f) + G<. (q, t - t')9(t' - i). 
Here G> , G< are 

Xj Xj 

+ e -i[^,x j (0)+iT Xj ](t-t') g ^ _^ 

G< j (c l ,t-t') = G> ] ( q ,t'-t), (7) 



where a; qiX3 (0) = ^q 2 + m\ x . + ReS Xj (q) + g]tp\, with S Xj . (q) the Xj field self-energy (recall 
that the field decay width T Xj is related to the imaginary part of the self-energy as T Xj (q) = 
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—lmH Xj (q)/(2u! cltXj )). Thus using Eq. (7) in Eq. (6), the explicit expression for Eq. (5) 
becomes 



d 3 q 



(x 2 ) ~ I — 

{Xj) 'J (2tt 



(2vr)3 2^.(0) 
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dt' 
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exp(-2r Xj 
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(2tt) 3 
sin[2w q , x .(0)|*-f , |] 



(8) 



For the second term on the RHS of Eq. (8), after integrating by parts with respect to t', it 
becomes 



f da U(t') - -4 ! 

J— oo L J 



d 3 q 
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exp(-2r Xi 
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(2tt) 3 

sin[2^ q , Xj (0)|t-t'|] = 



9) 



d 3 q 



1 



+ 9 



dt'<p(t')<p(t') 



(2vr) 3 4u. 
d 3 q 



q*(0) 



(27T) 



exp(-2r x3 |t-t / |) 



X 



{^ x /0)cos[2^ q , x /0)|t-t / |] + r x sin[2a; q , x /0)|t-f|]} 



(9) 



The first (local) terms on the RHS of both Eqs. (8) and (9), when pertubatively expanded 
in the coupling constant lead to quantum corrections from the Xj-fields to and A, to 
order and gj, respectively These corrections are divergent but are renormalized by the 
usual procedure of adding mass and coupling constant counter-terms. The second term on 
the RHS of Eq. (9) is responsible for dissipation. In this study, we are interested in the 
regime where ip(t) changes slowly relative to the relaxation time, in this case set by T Xj , 
which means the adiabatic approximation is valid. Under this approximation, similar to 
the treatment in [7], a Markovian, or equivalently time local, treatment can be used, which 
amounts to a derivative expansion of the field tp(t) and in which the leading ip term only 
is retained. After implementing this approximation and substituting Eq. (9) back into Eq. 
(8), we obtain 



(x 2 



d 3 q 



9 2 ^x 3 



(27r)32u qiX .(t) 



CO 



1,Xj 



(t) + n 



Xj 



(10) 



In the above, note we have conveniently reintroduced the time dependence back into the field 
frequencies and when they are perturbatively expanded to order gj, the above mentioned 
mass and coupling constant corrections are correctly reproduced. 



5 



An analogous expression to Eq. (10) also follows for (0 2 ). Note however that for an initial 
(at t — to) zero temperature bath and for fields $ and Xj satisfying the mass constraint 
m Xj > 2m^, fe > rritf,, there only will be decay channels for Xj into V'fe particles. As a result, 
it implies T^q) = and F^ k {q) = 0, while we have that 



ml I 4m? \ 3/2 



k=1 OTVUJq^ \ IIL X . J 



As such, in the adiabatic regime, dissipation will only involve the decay of Xj particles. 
The other two averages in the EOM, ((f) 3 ) and (4>Xj) can a l so be worked out in the linear 
response approach, and their leading contributions are at two-loop order [7]. Here, we will 
not consider them but restrict our calculation to leading one-loop order for simplicity. In this 
case, the only contribution to dissipation is Eq. (10), and this effect already will be adequate 
to demonstrate considerable radiation production from our model Lagrangian. Substituting 
Eq. (10) back into the effective EOM, Eq. (2), the second term on the RHS of Eq. (10) leads 
to a dissipative term in the EOM and the first term leads to $ mass and coupling constant 
divergent corrections, that can be renormalized as usual by the introductions of counterterms 
in Eq. (1). This renormalization procedure is standard and will not be further addressed. In 
our final expressions, all mass parameters, m 0( p, mo Xj ,mo^ k , and coupling constants, X,gj,hkj 
are then taken as the renormalized ones. The renormalized effective EOM for tp(t) that 
finally emerges can be written as 

0+^+^ = 0. (12) 

In the above equation, we have included in V e g the quantum renormalization corrections to 
the mass and coupling constant for the $ field, which are exactly the same as found in the 
calculation of a constant background (^-field effective potential. The dissipation coefficient 
7) (tp) in Eq. (12) comes from performing the momentum integral in Eq. (10) and using (11) 
to give 



-1/2 



( \ 2/ + \ V~* a X,W \ Xj 1 X,VJ / 10 \ 



3=i 32vr (2^ + ^ + 2m 



where = e£i h\ j:X m 2 x . (l - ^m\Jm\ J 3/2 /(8tt) and m Xj in Eq. (13) denote the field 
dependent masses, m x = m Xj (ip) = rn^ x . + g^ip 2 (t). The dissipative mechanism Eq. (13) 
overcomes an underlying impediment to realizing robust warm inflation in the finite tem- 
perature calculations [4,6], where all mass scales were constrained by the temperature. In 
sharp contrast, a key feature about the dissipative mechanism of this paper is that irrespec- 
tive of the magnitude of ip and m Xj , dissipation occurs unchecked by these severely limiting 
constraints. 

For the dissipative mechanism derived in this letter to be applicable to warm inflation, 
there must be some control in determining the quantum corrections in V e g in Eq. (12). This 
is required mainly since, similar to supercooled inflation, in the warm inflation case also, 
treatment of density perturbations requires an ultraflat potential [2,3,15]. However, there 
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are one-loop quantum corrections to the T = effective potential arising in the Lagrangian 
Eq. (1) from the self-interaction of the 0-field and from its interactions with the ^-fields, 
which give [16] 

1 r H 3 k ( Nx \ 



where E m<j> = ^/k 2 + + \(p 2 /2 and E mx . = ^/k 2 + m^. + g 2 (p 2 . To obtain the desired 

ultraflat potential, it requires A to be tiny with ~ A(^ 2 /2. In this regime, the contribution 
from the E m<t> term above is negligible. However, since in general we want gf ^> A, the one- 
loop contributions from the E m terms lead to corrections ~ gfip 4 in V e g and thus would ruin 
the flatness of the potential. Operationally these one-loop contributions can be controlled by 
adding to the Lagrangian Eq. (1) fermionic "partners" ip x to the x~ fields, with one ^x- field 
for every four x-fields and coupling only to the $-field as Ya=( A of^f^f^- The one-loop 
quantum corrections to the effective potential from these terms will yield [16] 

where E m ^ x = ^Jk 2 + (m ^x + gfip) 2 . In particular, this fermionic contribution has the fa- 
miliar opposite sign to the bosonic contribution. Thus with appropriately tuned parameters 
<?j, gf and with zero explicit masses m ^ x , = m 0Xi = 0, the one-loop quantum corrections to 
V e g cancel to all orders in gi, gf. This modification simply is mimicking super symmetry. For 
realistic model building, the mechanism derived in this letter must be examined in actual 
SUSY models, where the choice gf 3> A of coupling parameters could be obtained naturally, 
but that will not be pursued here. 



III. ALTERNATIVE DERIVATION OF DISSIPATION - OPERATOR 

FORMALISM 



For completeness, here an alternative derivation of dissipation is presented using the 
canonical approach and following the formalism developed in [12,14]. In this approach, the 
fields 0, x an d are expressed in terms of their mode decompositions and dynamics is 
determined with respect to the mode operators. Thus, for example for the Xjfat) field this 
means 

ttCM) = / ~ ^f q -771 [<W0e- iq " x + < {t)e*~] . (16) 
J (27r)2[2w q , Xj (t)]2 

Since there is a time dependent background field (p(t), this induces time dependence in the 
frequencies and so in the creation/annihilation operators of the <fi and Xj fields. In the 
analysis that follows, we will focus on the Xj fields, with similar considerations carrying over 
for the field. 
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The time dependent Xj - frequency in Eq. (16) is given by Wq iXj (t) = [q 2 + vn^ x . + 



£,V(*)] 1/2 . From Eq. (16) it follows that 



<*?> = / 



d 3 q 



2x q , x .(*) + 2Re[y q , x .(0] + l 



(17) 



where Xq iXj (t) = ( a q , Xj (£) a q,xj(^)) 1S the particle number density and Vq, Xj {t) = 
( a q,Xj W a -q,Xj * s the off-diagonal correlation. 

From the field equation for Xj an d Eq. (16) we can deduce the equations satisfied by 
x qiXj and i/q tXj - Taking also into account the possibility that the field Xj can decay into 
lighter fields with a decay rate T Xj (g) as already given in Eq. (11), %q, Xj and y q , Xj can be 



shown to satisfy the coupled differential equations [12,14] 



CO, 



x 



q.Xj 



y<i,xj 



co, 



Rey q , Xj , 



CO, 



co, 



q.xj 



iT xM) 



x q,Xj + 2 



- 2i 



^q.X, -^X 3 (^) 



Xj 



(18) 



A solution for Eq. (18) can be found in the quasi- adiabatic regime as follows. Let us consider 
the case of a slowly changing configuration tp(t). We can therefore suppose that the number 
of produced particles at time t is Xq jXj (t) 1. Consequently we also have that w qiXj and its 
time derivative slowly change. We then find for y qiXj in Eq. (18) the result 



S/q*,-(*) = 



which in the limit t ^> T^ 1 yield 



Uq, Xj {l -exp 


- 2 * (^q,X 3 -^X,)^ 


} 


4 


^q.X, ■ - iV X 3 (O) 


2 



(19) 



9) . 
2 W " 



Xj 



CO 



q.Xj 



r 2 

Xj 



(20) 



Using Eq. (20) in Eq. (17), once again we get Eq. (10), from which the effective 
EOM Eq. (12) follows. A shortcoming of this approach is that interactions are added to 
the set of Eqs. (18) in a somewhat ad-hoc way. This point was discussed recently in [8], 
where the complete kinetic equations where derived for the single field self-interacting c/> 4 
model. Nevertheless, the final answer from the approach of this section agrees with that 
from the Lagrangian based approach of the previous section, where interactions can be added 
consistently through the appropriate set of Schwinger- Dyson equations for the propagators 
[7]. Thus it suggests the results by this canonical approach are acceptable, but missing gaps 
in the formalism of [12] must still be resolved. For our purposes, due to the importance of 
the dissipative mechanism studied in this letter, we felt it was important to point out the 
agreement between independently developed formalisms, even if there remain shortcomings 
in one of them. The practical significance of the results in this letter provide motivation to 
address these difficult problems in the course of future work. 
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IV. PHYSICAL INTERPRETATION AND AN EXPLICIT APPLICATION 



We now turn to an application of the equations derived above, using an explicit set of 
model parameter values, which are consistent with simple inflationary models. But before 
that, let us address briefly the physical interpretation of dissipation in Eq. (12). 

We note that the evolving background field ip(t) changes the masses of the Xj bosons. As 
a consequence, the positive and negative frequency components of the Xj-fields mix. This in 
turn results in the coherent production of Xj particles which then decohere through decay 
into lighter -?/> fc -fermions. This picture can be confirmed by checking energy balance. This 
is done by examining the time evolution of the x_,-particle number density. For this, their 
number density is expressed in terms of time dependent creation and annihilation operators 
as Af = J2j( a l Cj {t) a xj('t)) ■ By relating the time dependent operators a> Xj (t) and a Xj (t) to 
the initial, time independent, creation and annihilation operators through a Bogoliubov 
transformation, the total particle production rate then can be computed in general. Thus, 
the time evolution of the total production rate is 

which using Eqs. (18) and (20), leads to 

Nx r rl^k a 4 Y 

_ajt g, 1% (22) 

It can now be checked from Eqs. (2), (10) and (12), that the above result, Eq. (22) is 
precisely equal to the vacuum energy loss rate, r]{p 2 , as obtained from the effective EOM, 
Eq. (12). 

Let us now examine the application of the results in this letter to warm inflation and also 
understand their significance. The scope of the present calculation is limited since dissipation 
at zero temperature necessarily implies a nonequilibrium state, which is evolving to some 
statistical state containing particles. Thus the estimates made below only give some idea 
of the magnitude of particle production. However, provided the magnitude is significant, 
as will be shown, it reveals that on scales relevant to inflation, quantum field theory with 
generic interactions has robust tendency to dissipate. For our estimates, we have set same 
all $ — x couplings g Xj = g as well all x ~ i> couplings, hkj = h. 

We are interested in overdamped motion for the inflaton tp(t), which requires i). = 
m 4,( i P) = m 0(/> + A<£> 2 /2 < ? ? 2 (v 9 )/4- The adiabatic approximation Eq. (3) requires ii). 
m^(ip) / r]((p) < T x . Although our derivation was for Minkowski spacetime, provided the 
time scale of microscopic dynamics is faster than the Hubble time scale, then within sub- 
Hubble length scales, this Minkowski spacetime calculation should be valid. For this to 



hold, it requires iii). H = ^/87rV r cff /3mp 1 ~ ^/87r(A/4!) Lp 4 /3m^ < T x , where m pl is the 
Planck mass. Also, so that the macroscopic motion of tp is governed by the dissipative term 
it requires iv). r](ip) > 3H. Thus combining all four of the above consistency conditions 
leads to parametric constraints. To obtain these, we will treat m x > a x ^, where from 
below Eq. (13) we have, by setting m x ~ g 2 <p 2 , o? ^ « g 2 h 2 N^(p 2 /(8ir), which thus requires 
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h 2 N^/(8n) < 1. In this regime, we have from Eq. (13) rj xs g 3 h 2 N x N^(p/(5127r 2 ) and 
from Eq. (11) T x « gh 2 N^p/(8iT). The parametric constraints that follow from the four 
conditions given above are respectively 



i. A < 

ii. A < 

iii. A < 



2(512tt 2 ) 2 

9 4 h A N x N 2 

2048vr 3 
9g 2 h*N 2 m\ x 



iv- A < • (23) 



64tT 3 (/3 2 

gWN 2 N 2 m 2 pl 
7r(512vr 2 ) 2 v? 2 

To yield large dissipation, we are usually interested in the regime where the couplings g, h 
are big. To remain within a well defined perturbative region, we will then further require 
that g 4 N x < 1 and h 2 N^ < 1 and will base our estimates on the upper bounds here. Also, in 
general ip ^ m p i, but to obtain the tightest constraints on A in (iii) and (iv), we will set this 
at the equality point. Under these conditions, we find for the constraints (i) - (iv) in Eq. 
(23) respectively A < min(10 -8 g 2 iV x , 10~ 5 , 10~ 3 g 2 , 10~ 8 g 2 iV x ). Recalling that constraints 
imposed by density fluctuations give typically A < 10~ 14 [2,3,15], we see that the above 
constraints introduce no stricter limitations. 

As shown in Eq. (22), radiation production is determined by 



dv es ... ,. ,._^ m2 M 



p r (t) = V (p)p 2 = — » V«(tp)-*^. (24) 
dp rj 

The zero temperature calculation should be valid for a time period ~ l/r x , in which time 
the magnitude of radiation produced is 

Pr (l/r x ) » V eS (<p)ml(y>)/(r)r x ) < K ff (^). (25) 

Based on Eqs. (11) and (13) and the above constraints on A, there is considerable freedom 
in choosing the ratio TZ = m 2 ./(^r x ) appearing in Eq. (25). Considering then an ultraflat 
potential, as necessary for observationally consistent density perturbations, which requires 
typical values of A ^ 10 -14 , this implies TZ ^$ 10 _10 /((7 4 /i 4 A^A^ x ). For unexceptional values 
of the perturbative coupling parameters, say g ~ h ~ 0.1, and small number of x an d ip 
fields, N x , N^p ~ 1 — 10, this leads to 1Z ~ 10 - *- 2-5 - 1 . Also note these parameters choices 
are consistent with the conditions on A given above Eq. (24). Thus for a typical scale for 
inflation, where the potential energy is at the GUT scale, Ves(p) 1 ^ ~ 10 15_16 GeV, it implies 
a generated radiation component which, if expressed in terms of temperature, is at the scale 
T ~ 10 13 ~ 16 GeV, and this is non-negligible. This is a significant result not only because the 
magnitude of produced radiation is large, but also because it emerges from a very generic 
interaction, scalar — > heavy scalar — > light fermions, which is very common in many particle 
physics models. Moreover, we expect similar robust radiation production for decay of the 
heavy scalars into gauge bosons. Finally, although we did this zero temperature calculation 
first simply due to its tractability, an interesting fact emerges for inflationary cosmology, 
that even if the initial state of the universe before inflation is at zero temperature, the 
dynamics itself could bootstrap the universe to a higher temperature during inflation. 
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V. EXTENSION TO EXPANDING SPACE-TIME 



The extension of this calculation is formally straightforward to Friedmann-Robertson- 
Walker (FRW) spacetime, ds 2 = dt 2 — a 2 (t)<ix 2 , where a(t) is the cosmic scale factor and 
t is cosmic time. In this case, the extension of Eq. (1), for the Lagrangian density of the 
matter fields coupled to the gravitational field tensor g^ u , is given by 



C = 



Ny r 



-9 -g^dfidu® 



2 4! 



-m 2 

2 



+E 

iVy, 

+E 

k=i 



m, 



4! 



Ny 



(26) 



where R is the curvature scalar and £ is the dimensionless parameter describing the coupling 
of the matter fields to the gravitational background. In the last terms involving the fermion 
fields, the 7 M matrices are related to the vierbein (where = e a e\ J r] a \ ) , with i] a b the usual 
Minkowskii metric tensor) by 7 M (x) = 7 a e^(x) [17], where 7 a are the usual Dirac matrices 
and = -{i/4)<j ab e v a V^e bv , with a ab = i/2[ 7 a , 7 6 ]. 

It is easy to show that the Lagrangian Eq. (26) in conformal time, t c , where dt = adt c , 
remains unchanged from Eq.(l) except that all masses obtain time dependence related to 
a(t c ) (see for example [14] for a similar problem). In particular, for the bosonic fields we 



have that m^(t c ) = ml Xj a 2 (t c ) — d 2 a/2adt 2 + £a 2 i?/2 and similar for the field, and for 
the fermionic fields m^ k (t c ) = m ^ k a(t c ). These time dependent parameters can be treated 
within the linear response formalism used in this letter. Moreover, since the time dependence 
is associated with a(t c ), it is easy to show that provided H < T x , the time dependence of 
the mass terms is slow relative to microscopic dynamics and thus an appropriate adiabatic 
approximation should be applicable. 

The observations made above are adequate to establish that, for the mechanism of central 
interest in this letter, the robust dissipative properties found above for Minkowski spacetime 
also will hold for expanding spacetime. However, the exact form of the effective <p-EOM 
is a more involved matter. The problem is there are three relevant time scales H, T Xj 
and y?/V) where for the slow-roll motion of interest, we seek solutions with ip/<p < H. 
Moreover, ultimately we require the evolution equation in cosmic time, and the relation 
between that and conformal time is in general very nonlinear. For example, for the case 
of prime interest, de Sitter space, t oc ln(l — bt c ). Thus power law ambiguities can have 
nontrivial relevance in relating between conformal and cosmic time, and such ambiguities 
are prevalent in adiabatic approximations and derivative expansions. This is a serious matter 
and to learn more about this mechanism in expanding spacetime beyond what already has 
been understood from the above Minkowski spacetime calculation requires application of 
more complete nonequilibrium methods, such as [18]. We will consider the details of this 
derivation in the FRW spacetime in a future work. 
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VI. CONCLUSIONS 



The relevance of the analysis in this letter extends beyond warm inflation, since the 
interactions studied here are exactly the same as found in supercooled inflation models. In 
fact, in the context of the model studied here, with couplings around the ones studied in the 
example of Sec. IV, reheating becomes irrelevant, since our analysis showed the model is 
inconsistent with supercooling in the first stage, and the entire dynamics is warm throughout. 
Thus, as originally suggested [2,15], warm inflation dynamics is inherently intertwined with 
the general problem of inflationary dynamics. 

Since the first principle results in this paper give support to the warm inflation picture, 
it is worth recalling here other features that also have made this picture compelling. First, 
warm inflation overcomes a conceptual barrier that the supercooled picture has never shaken 
away, which is that in warm inflation there is no quantum-classical transition problem, since 
the macroscopic dynamics of the background field and fluctuations [15] are classical from 
the onset. Second, in warm inflation models, in regimes relevant to observation, the mass of 
the inflaton field is typically much larger than the Hubble scale, thus these models do not 
suffer from what is sometimes called the "eta problem" . Finally, accounting for dissipative 
effects may be important in alleviating the initial condition problem of inflation [19,20]. 

The emerging picture is that warm inflation remains a hopeful direction toward a com- 
plete and consistent dynamical description of the early universe. However, considerable work 
remains in understanding the quantum field theory of this picture. Two areas were already 
identified in the paper. One is resolving the gaps in the canonical dissipative formalism of 
[12], thus permitting this approach to be a viable cross-check to the Lagrangian approach. 
The other area is a full investigation of the dissipative formalism in expanding spacetime. 
Beyond this, the more difficult problem is extending the adiabatic contraints in the present 
formalisms to treat nonequilibrium conditions. Steps along this direction already have be- 
gun, using operator methods [9] and the even more ambitious attempt in [8] to derive the 
Boltzmann-like kinetic equation for interacting fields. 
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